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Capitulo 9. Fotometria astronòmica 

Estos apuntes para el curso AS750 

se encuentran en 

http;//www.das.uchile.cl/~mhamuy 

/courses/AS750/AS750.html 

Leetura Obligatoria: 

Lectura Sugerida: Gutiérrez & 
Moreno (Capitulos III y VII) 

Sección 1. Introducción 

Idealmente uno quisiera conocer la 
distribución espectral de las 
fuentes astronómicas, pero la 
dispersión de la luz resulta cara 
en tiempo de telescopio. Para 
obviar este problema el observador 
selecciona un amplio rango 
espectral (banda de paso) para lo 
cual utiliza un filtro con lo cual 
obtiene un promedio del flujo en 
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dicho rango. Este tipo de medición 
se denomina fotometria 
astronòmica, cuyo objetivo 
consiste en medir el flujo de las 
fuentes a través de bandas de paso 
estàndares pre-establecidas. 

Recordemos del Capitulo 2 que la 
magnitud de banda ancha se define 
corno: 


mj = 


- 2 . 



n° k Sx dX+Cj = 



X 

he 


Sx dX 


en que: 

on.° es el flujo de fotones (# s"' 

A 

À 1 cm" 2 ) fuera de la atmosfera 
terrestre 

of 0 es el flujo de energia de la 

fuente (erg s 1 À" 1 cm’ 2 ) fuera 
de la atmosfera terrestre 
oS. es la función de transmisión 

A 
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que define la banda de paso 
estàndar 

oC. es una constante aditiva 

A 

arbitraria para cada banda de 
paso 

Como el observador ve un flujo 
modificado por los efectos de la 
extinción atmosfèrica y por su 
propia respuesta instrumentai, su 
medición resulta en una magnitud 
instrumentai : 




/a Yc E ? A dX 


en que E. x es la transparencia de 

A 

la atmosfera (que es función de À 
y de la masa de aire X), S' es la 

banda de paso instrumentai y A es 
el area colectora del telescopio. 

El desafio del observador consiste 
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en remover los efectos 
instrumentales y de extinción y 
convertir sus magnitudes 
instrumentales (o naturales) en 
magnitudes en el sistema estàndar 
(magnitudes estàndares). 

Sección 2. Extinción atmosfèrica 

La luz de las fuentes astronomicas 
se atenua al cruzar la atmosfera 
debido a absorción y scattering 
por lo cual se hace necesario 

corregir este efecto. Si f.° (erg 

A 

s" 1 À" 1 cm" 2 ) es el flujo de una 
fuente fuera de la atmosfera el 
flujo recibido en la superficie de 

la Tierra (f.) es menor que f.°. 

Para ver còrno se atenua el flujo 
partamos por la ecuación de 
transferencia radiativa: 
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dfx = -fx kx p(h) ds 


en que k. (cm 2 g 1 ) es la opacidad 

A 

de la atmosfera, p (g cm" 3 ) es la 
densldad de la atmosfera a una 
altura h sobre la superficie, y ds 
es el elemento de distancia entre 
la estrella y el observador. 

Suponiendo una atmosfera de capas 
planas y un objeto celeste a una 
distancia cenital z: 


cenit 



se puede ver que : 
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ds = — sec z dh 


con lo cual: 


- 4 — = k\ p(h ) sec 2 dh 

Jx 


La solución es: 



sec z 


J 0 H p(h) dh 


en que H es la altura total de la 
atmosfera medida desde el 
observatorio. Si defino corno la 

masa de aire por unidad de àrea en 

la vertical (g cm" 2 ) y p corno la 
masa de aire por unidad de àrea en 
la dirección de observación, 

r H 

li {) = / p( h) dh p = pò sec z 

Jo 


y si defino la masa de aire X (sin 
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dimensiones) en la dirección de 
observación. 


X = — = sec z 
Vo 


la ecuación de transferencia 
queda: 


fx = fi e"** mX 

El flujo decae exponencialmente y 
el decaimiento crece con la masa 
de aire. 

La masa de aire es 1 en la 

vertical (z=0), y 2 para z=60°. La 
aproximación de capas planas es 
valida cerca del cenit. Tornando en 
cuenta la curvatura de la 
atmosfera la masa de aire es: 


X = sec z - 0.00181G7(sec z- 1) - 0.0028750(sec 2 - l) 2 ... 
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Podemos definir la transparencia 
E. de la atmosfera en la vertical 

A 

(X=l): 


= 1 ) 


I 




(Ex < 1 ) 


con lo cual vemos que la 
transparencia depende de la 
opacidad 



g-*'À M) 


El flujo de la estrella para una 
masa de aire X es: 

h = fi Et 


Sección 3. Respuesta instrumentai 
La banda de paso instrumental, S 1 . 

A 

es el resultado de la atenuación 
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de fotones por todos los elementos 
ópticos del sistema de detección. 
Como minimo se deben considerar la 
reflectividad del espejo(s) (M. ), 

la transmisión del filtro(s) (F.) 

A 

y la eficiencia cuàntica del 
detector (Q.). Si consideramos un 

A 

telescopio tipo Cassegrain habrà 
dos reflexiones, por lo que la 
banda de paso sera: 

ST X = Mi Q x Fx... 

En generai habria que agregar a 
està ecuación factores adicionales 
que den cuenta de todos los 
elementos ópticos a lo largo del 
camino óptico hacia el detector. 
Veamos algunos de estos elementos: 

oEficiencia cuàntica del 
detector: 
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El primer detector astronòmico 
fue el ojo humano. Desde fines 
del siglo XIX las placas 
fotograficos comenzaron a 
utilizarse en astronomia y 
produjeron una revolución porque 
permitian tener un registro 
permanente de las observaciónes. 
Las emulsiones tenian 
eficiencias cuànticas bajas, 
entre 1% y 5%. Alrededor de 1940 
se comenzaron a usar los tubos 
fotomultiplicadores con 
eficiencia de hasta 90% aunque 
con un solo canal. A comienzos 
de 1980 se comenzaron a usar los 
detectores bi-dimensionales en 
longitudes de onda ópticas 
(CCDs) y unos pocos anos mas 
tarde los primeros arreglos para 
la detección de luz en el IR. 
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ìMtvdength [A] 


oReflectividad de espejos: 

Habitualmente los espejos se 
recubren con delgadas capas de 
aluminio aunque hay nuevos 
materiales mas reflectantes 
dependiendo de la longitud de 
onda de interés. 
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wavelength [A] 


oFiltros: 

Los filtros se fabrican con 
distintos materiales ópticos con 
el fin de aislar ciertas zonas 
espectrales de interés. Cuando 
el ancho de la banda a media 
altura es mayor que 300 À se 
habla de fotometria de banda 
ancha, entre 100-300 À es 
fotometria de banda intermedia, 
y para anchos menores a 100 À se 
trata de fotometria de banda 
angosta. 
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A continuación se muestra uno de 
los sistemas fotométricos mas 
usados: el sistema de banda 
ancha UBVRIJHKLM establecido por 
Johnson entre 1953-1959. 



el espectro de Vega 

Acà se muestra el sistema de 
Sloan introducido en 1996 que 
està siendo muy empleado. Es un 
sistema de banda ancha hecho de 
filtros interferenciales los 
cuales permiten definir muy bien 


13 of37 


05/07/2010 04:17 PM 






























Clase 9 


http : //www. das. uchile. cl/~mhamuy/courses/AS 750/tema9... 


los limites de cada banda 
(excepto el filtro z'), y una 
mejor reproducibilidad. 



Filtros del sistema de Sloan y 
el espectro de Vega 

Los filtros se caracterizan por 
la longitud de onda equivalente 
y/o longitud de onda efectiva : 

reo reo 

/ S\ A d\ = X eq / S\ d\ 

J[) J{) 
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rCO rCO 

/ fx S\ A d\ = A,,/; / f\ S\ 

■In Jn 


d\ 


Caracteristicas de los filtros 
de Johnson y de Sloan: 


FILTRO 

A 

eq 

FWHM(À) 

U 

3652 

524 

B 

4468 

1008 

V 

5505 

826 

R 

6581 

1576 

I 

8059 

1543 

j 

12400 

2300 

H 

16500 

3100 

K 

21600 

3300 

L 

35000 

4000 

M 

47500 

4600 

u 1 

3522 

634 

g' 

4803 

1409 

r 1 

6254 

1388 

i 1 

7668 

1535 

z 1 

9114 

1409 


El sistema UBVRI està definido 
en Bessell 1990 (PASP, 102, 

1181) 

Los filtros JHK estàn definidos 
en Persson et al. 1998 (AJ, 116, 
2475) 
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Los filtros L y M estàn 
definidos en Bessell and Brett 
1988 (PASP, 100, 1134) 

El sistema de Sloan està 
definido en Fukugita et al. 1996 
(AJ, 111, 1748) 

oTransmisión de otros elementos 
ópticos en el camino óptico: 
ventana del CCD, correotor de 
campo, etc. 

Los datos para producir estos 
gràficos estan ACA . 

Sección 4. Corrección por 
extinción atmosfèrica 

La magnitud instrumentai que el 
observador mide es: 



y lo que queremos obtener es la 


16 of37 


05/07/2010 04:17 PM 



Clase 9 


http : //www. das. uchile. cl/~mhamuy/courses/AS 750/tema9... 


magnitud fuera de la atmosfera 
(E À =1), i.e., 



7 A Si dX 
he A 


La magnitud instrumentai se 
relaciona con la magnitud fuera de 
la atmosfera por: 

rn^ = tra® + Kj X 


en que K. es el coeficiente de 

J A 

extinción [en unidades de magnitud 

x (masa de aire) 1 ] el cual 
corresponde a la transparencia 
media para la banda de paso: 


Kj = -2.5 log Ej = -2.5 log 


r jr n ^ % dx \ 
L r n a si dx j 


l/.V 
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obtendremos la recta de Bouger: 



Para determinar el coeficiente de 
extinción K. debemos restringir el 

A 

rango a X=l-2 ya que a masas de 
aire mayores la aproximación 
lineai ya no es vàlida. Està 
medición sólo se puede hacer 
cuando la extinción se mantiene 
estable en el tiempo y en 
distintas direcciones durante la 
noche, en cuyo caso se habla de 
"noche fotomètrica". Para que se 
cumpla està condición al menos la 
noche debe estar totalmente 
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despej ada. 

Coeficientes de extinción medidos 
en Las Campanas en el sistema de 
Johnson y Sloan: 


FILTRO 

A 

K\ [mag (masa de aire) ] 

U 

0.57 

B 

0.29 

V 

0.13 

R 

0.10 

1 

0.08 

J 

0.11 

H 

0.05 

K 

0.08 

u' 

0.51 

g' 

0.20 

r' 

0.11 

i' 

0.06 


Notar que la extinción disminuye 
fuertemente con la longitud de 
onda. 

Sección 5. Transformación de 
magnitudes instrumentales al 
sistema estàndar 

Una vez determinada la extinción 
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podemos calcular magnitudes 
instrumentales fuera de la 
atmosfera. Bien, pero las 
magnitudes instrumentales varian 
de instrumento a instrumento. Para 
poder comparar magnitudes 
obtenidas con instrumentos 
diferentes se define un sistema 
estàndar al cual deben convertirse 
las magnitudes instrumentales. 

Si S. es la banda de paso del 

A 

sistema estàndar, la magnitud en 
el sistema estàndar (o magnitud 
estàndar) se define corno: 


™A = 


-2.5 



«S& dX+Cj = 



X 

he 


S\ dX 


en que C. es el punto cero de la 

A 

escala de magnitudes y difiere 
para cada banda de paso. 


Notar que la banda de paso 
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instrumentai no es igual a la 
banda estàndar (i.e. S' * S,), 

A A 

por lo cual la diferencia entre la 
magnitud instrumentai fuera de la 
atmosfera y la magnitud estàndar 
es : 


mj — m? = —2.5 log 


j,r f" x d\ 

f~n A A S{ dX 



El término con el àrea es 
independiente de la estrella 
medida y se puede combinar con C. 

A 

en una constante aditiva ZP. a ser 

A 

determinada observacionalmente: 


rnj— nij = —2.5 log 


Io°° fx A dX 

r fx A SI dX 



El término con la razón de los 
flujos, sin embargo, depende de la 
distribución espectral de cada 
estrella: 
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La ecuación de tranformación al 
sistema estàndar se puede 
escribir: 


mj = rrij + (C r l\ x color) + ZPj 


en que CT es el término de color. 
Para determinar CT. y ZP. debemos 

A A 

observar estrellas estàndares de 
distintos colores. Como hay dos 
incógnitas al menos necesitamos 2 
estrellas. Lo ideal es tener 20 
estrellas y usar un ajuste de 
cuadrados minimos. 


Si S \ = 
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valor medido de CT. es un 

A 

diagnòstico de cuàn parecido es el 
sistema instrumentai al sistema 
estàndar (mientras mas chico 
mejor) . 

Ejempio de ecuaciones de 
transformación: 

U = u — K V X + CT v (u - b) + ZPu 
B = b- I< B X + CTg(b- v) + ZP D 
V = v - K V X + CT v (b -v) + ZP V 
R = r — I< R X + CTr(v -r) + ZP R 
I = i- I<iX + CT[(v - i) + ZPi 

Catàlogos de estrellas estandares 
en el sistema UBVRI: Landolt 1992 
(AJ, 104, 340) 

Catàlogos de estrellas estandares 
en el sistema JHK: Persson et al. 
1998 (AJ, 116, 2475) 

Catàlogos de estrellas estàndares 
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en el sistema Sloan: Smith et al. 
2002 (AJ, 123, 2121) 

Sección 6. Ejempio de fotometria 
CCD en el sistema UBVRI 

Noche del 30/09/2002 con el 
telescopio de 1 m de Las Campanas. 
Se observaron 5 campos de 
estrellas estàndares de Landolt y 
varios campos de objetos de 
programa. Se usò el paquete 
"apphot" en IRAF para la reducción 
de los datos. Està figura muestra 
el campo T Phe y las cuatro 
estrellas estàndares: 
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* 


B 

* 

* 


P 


«Medición de magnitudes 
instrumentales: 


Se usò la tarea "phot" para 
medir los fotones (electrones) 
registrados en el CCD de las 
estrellas de Landolt usando una 
apertura de radio de 7". Està es 
la misma apertura usada por 
Landolt en el establecimiento 
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del sistema fotomètrico UBVRI. 
Para restar el cielo se 
seleccionó un anillo de 2" de 
ancho a una distancia de 8" de 
cada estrella. Acà se muestra el 
perfil radiai de una estrella y 
los pixeles en el anulo del 
cielo: 
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La tarea "phot" calcula el cielo 
medio en el aniIlo, resta este 
cielo de cada pixel contenido en 
la apertura y suma la 
contribución de los fotones 
registrados en cada pixel dentro 
de la apertura. La tarea se 
encarga de normalizar el nùmero 
de fotones por el tiempo de 
exposición y de calcular una 
magnitud instrumentai. 

oErrores instrumentales: 
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Supongamos una apertura con n # 
pixeles, un ànulo con n 

s 

pixeles. Si p son los fotones 

registrados en el pixel i debido 

a la estrella y s. los fotones 

1 

registrados en el mismo pixel 
debido al cielo (i:l ->• n # ) , los 

errores en cada cantidad son 

Pi± Vpì 

Sj i y/si ~ 8 ± Va 

El cielo medio se obtiene de los 
j pixeles del anulo (j:l ->• n ) 


S = 


E, 


± 


n, 


vi 



n s 


Como no podemos separar los 
fotones de la estrella y del 
cielo en el pixel i, el flujo 
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neto de la estrella en el pixel 
i lo obtenemos restando el cielo 
medio medido en el ànulo: 

fi = Pi + s i — s 

y el flujo neto total de la 
estrella es: 


/ = 7>. + Si — «) 

i 


El error en f se obtiene sumando 
en cuadratura: 

a ì = + a l + a *+ ron2 ) = ip ‘ + 5 '+~+ r °n 2 ) 

i i S 

en que ron es el ruido de 
lectura del CCD. En el caso n 

S 

» 1 : 


o\ = / + n*s(l H-) + n*ron 2 

J ' n/ 
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El error relativo en f es: 



1 n*(s + ron 2 ) 


s + r- 


Notar que a /f decrece cuando f 

crece. Dependencia de a /f con 

: el primer término crece 

cuando disminuye porque se 

pierde flujo al achicar la 
apertura, mientras que el 
segundo término decrece cuando 
disminuye. Hay un punto 

óptimo en que a /f es minimo y 

este punto depende de la 
distribución radiai de f, i.e., 
el seeing. 

La relación entre a,/f y a se 


m 
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obtiene usando la fòrmula de 
propagación de errores: 


m = —2,5 log / + eie 


Sm 


2,5 


Sf f In 10 


2,5 Uf ai- 

= sr tòt ^ Lm 7 

Notar que si a =0.01 mag ->■ 

J m 

(a f /f) =1%, -► (S/N) =100 

oTransformación al sistema 
estandar: 

La tarea "fitpar" usa las 
magnitudes instrumentales y las 
magnitudes estàndares de Landolt 
para resolver el coeficiente de 
extinción, término de color y 
punto cero de la transformación 
usando un ajuste de cuadrados 
minimos pesados. 
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V = v - K v x + CT v (b - v) + ZPv 

Estos gràficos muestran los 
residuos de la transformación 
asumiendo que el coeficiente de 
extinción es cero y luego 
asumiendo que el término de 
color es cero: 


NOAO/IRAF V£ +12 + 2-EXPORT mhamuyetennis Tue 17:30:07 ll-Apr-200É 
low_rej=3 + , high_rej=3+, nreject=3^ grow=0 + 
tota1=25, rejeoted=0, deleted=3, RMS=0 + 02801 
tolerance=3+000E-5, maxiter=15, iterations=ll 
V = v-kv^xv+ctvM(v-i>+zpv 
Solution converged 
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NOAO/IRAF V2 +12 + 2-EXPORT mhamuyOtennis Tue 17:33:39 ll-Apr-200É 
low_rej=3 + , high_rej=3+, nreject=3, grow=0 + 
tota1=25, rejected=Q, deleted=3, RHS=INDEF 
tolerance=3+00GE-5, maxiter=15, iterations=15 
V = v-kv^xv+ctvM(v-i>+zpv 
Solution did not converge 



Ahora se muestran los residuos 
de la transformacion ajustando 
el coeficiente de extinción y el 
término de color: 


NOAO/IRAF V2 +12 + 2-EXPORT mhamuyOtennis Tue 17:35:15 ll-Apr-2GGÉ 
low_rej=3 + , high_rej=3 + , nreject=3, grow=G + 
total=25, rejected=0, deleted=3, RHS=0 + 01105 
tolerance=3+000E-5, maxiter=15, iterations=4 
V = v-kvMxv+ctvM(v-i>+zpv 
Solution converged 
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HOAO/IRAF V£ +12 + 2-EXPORT mhamuyetennis Tue 17:34:10 ll-Apr-200É 
low_rej=3 + , high_rej=3 + , nreject=3, grow=0 + 
total=25 r rejected=0^ deleted=3.r RHS=0 + 01105 
tolerance=3+000E-5, maxiter=15, iterations=4 
V = v-kvHxv+ctvH(v-i^+zpv 
Solution converged 



Este ajuste arroja una dipersión 
(desviacion estàndar) de 0.011 
mag. Cuando la dispersión 
sobrepasa 0.03 mag se sospecha 
de la calidad fotomètrica de la 
noche. 

oMedición de magnitudes para 
objetos de programa: 

Los objetos de programa suelen 
ser débiles por lo cual no 
conviene medir sus flujos a 
través de aperturas muy grandes. 
Con el fin de minimizar el 
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ruido, se mide el objeto de 
interés con una apertura pequena 
y se aplica una corrección de 
apertura obtenida a partir de 
algunas estrellas brillantes del 
campo (corno las indicadas 
abajo) . 



35 of37 


05/07/2010 04:17 PM 




Clase 9 


http : //www. das. uchile. cl/~mhamuy/courses/AS 750/tema9... 


determinar la curva de 
crecimiento, definida corno la 
diferencia en magnitud entre 
aperturas concéntricas, en 
función del radio de la 
apertura. Està curva mide cuando 
flujo se pierde al achicar la 
apertura. 



Finalmente, la magnitud 
instrumentai del objeto de 
interés se convierte al sistema 
estàndar usando la ecuación 
anterior. 
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